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We study the effect of quantum statistics on the arrival time distribution of quantum particles 
computed through the probability current density. It is shown that symmetrization or asymmetriza- 
tion of the wave function affects the arrival time distribution for even freely propagating particles. 

In order to investigate the effect of statistics on the weak equivalence principle in quantum mechan¬ 
ics (WEQ), we then compute the mean arrival time for wave packets in free fall. The violation of 
WEQ through the effect of statistics on the mass-dependence of the mean arrival time is clearly 
exhibited. We finally evaluate the effect of spin on the violation of WEQ using a different approach 
by including an explicit spin-dependence in the probability current distribution, and compare it 
with the approach using particle statistics. Our results show WEQ re-emerges smoothly in the limit 
of large mass. 
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I. INTRODUCTION 

The statistics of elementary particles is a key underlying ingredient of quantum field theory which has been most 
successful in explaining three, viz. electromagnetic, weak and strong nuclear forces, of the four fundamental forces in 
nature. On the other hand, the explicit role of quantum statistics in physical phenomena such as Bose Condensation 
[l| is also being increasingly displayed in recent years. Practical applications of such statistics is potentially wide, 
ranging from optical lattices 0 to quantum information Q. In most practical manifestations of the role of Bose 
statistics, a crucial role is played by the measurement of the time of flight distribution Q. For example, the coherence 
properties of a condensate of ultracold atoms is imprinted on the time of flight distribution when the trapping potential 
is released and the particles undergo free fall. Thus, the effect of gravity on the motion of quantum particles seems 
to play an important role in the inference of attributes related to the statistics obeyed by them. 

The motion of freely falling particles is traditionally taken to conform to the weak equivalence principle (WEP) of 
gravitation which states that all sufficiently small objects fall with the same acceleration independent of their mass 
and constituent in a uniform gravitational field. WEP is regarded to be a fundamentally classical and local principle. 
Study of the equivalence principle in quantum mechanics has evoked enough interest. A statement of the principle in 
quantum mechanics is as follows: ’’The results of experiments in an external potential comprising just a sufficiently 
weak, homogeneous gravitational held, as determined by the wavefunction, are independent of the mass of the system” 
0. This assertion is also called the weak equivalence principle of quantum mechanics (WEQ). 

Different approaches have been used to study the possibility of violation of weak equivalence principle in quantum 
mechanics, starting from the prediction of mass-dependence of the radii, frequencies and binding energy of a particle in 
an external gravitational field @ . Experimental probes into possibility of violation of the weak equivalence principle in 
quantum mechanics have been attempted using interference phenomenon associated with the gravitational potential 
in neutron and atomic interferometry experiments 0, @|- However, the phase shifts observed in such experiments 
0 may be expressed in terms of directly observable quantities which turn out to be independent of mass thus 
ruling out the violation in this context. Moreover, atomic interferometry experiments @ need to be performed with 
different atomic species in order for a rigorous experimental test of the weak equivalence principle at the quantum 
level. A gedanken experiment studying the free fall of quantum test particles in a uniform gravitational field predicts 
mass-dependence of the time of flight distribution |10l |. Another approach using a model quantum clock predicts 
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violation of WEQ in the vicinity of the turning point of classical trajectories [TT|. An explicit mass dependence of 
the position probabilities has been shown for quantum particles projected upwards against gravity around both the 
classical turning point and the point of initial projection using Gaussian [l^ and non-Gaussian [l^ wavepackets. 
Moreover, explicit mass dependence of the mean arrival time at an arbitrary detector location has also been predicted 
for a Gaussian [l^ wave-packet under free fall. Such mass-dependence may be enhanced using suitably chosen non- 
Gaussian wavepackets [l3|. Thus, both the position probabilities and the mean arrival time show the violation of 
WEQ for quantum particles described by Gaussian and non-Gaussian wavepackets. It has been discussed violation of 
WEQ arises as a consequence of the spread of wave packets, the magnitude of which itself depends on the mass. An 
illustration of this effect has been presented based on Bohmian trajectories [l3l| . 

The above studies have established the violation of WEQ in single particle quantum mechanics. However, the effect 
of quantum statistics on the violation of WEQ has remained uninvestigated till date. As pointed out earlier, several 
important phenomena based on quantum statistics is experimentally revealed through the measurement of time of 
flight of quantum particles in free fall. Consideration of quantum mechanical effects on such time of flight distributions 
beyond the standard semi-classical analysis could reveal interesting observational effects, as has been debated earlier 
in the literature H- In the present work we are interested to examine the effect of quantum statistics on the WEQ. 
With the above motivation in this paper we study the effect of statistics on the arrival time distribution of a system 
of freely falling wavepackets consisting of two identical particles. 

Our analysis is based on the probability current approach for computing the mean arrival time distribution of 
wavepackets [l^. Before embarking on the study of WEQ, we will first consider the free evolution of a system of 
two identical particles to investigate if there is any effect of Bose-Einstein (BE) or Fermi-Dirac (ED) statistics on the 
arrival time distributions. Quantum statistics is fundamentally related to the spin of the quantum particles. It has 
been shown earlier that inclusion of spin enables a unique determination of the relativistic probability current [l^ [T^ 
leading to an explicit spin-dependent effect at the non-relativistic level for freely evolving fermions [l^ as well as 
bosons [l^ . Here we present a calculation of the effect of such explicit spin-dependence on the violation of WEQ for 
a system of two identical particles in free fall. 

The plan of this paper is as follows. In the next section we present the framework for studying the arrival time 
distribution for a system of two particles using the probability current approach. We show that symmetrization or 
anti-symmetrization of the wavefunction leads to differences from the case of Maxwell-Boltzmann statistics when there 
is overlap of the wavepackets. In section HI the effect of BE and ED statistics on freely propagating wavepackets is 
studied. In section IV we address the main motivation of the paper by studying the violation of WEQ separately 
using two different schemes, viz, (a) effect of statistics on the arrival time distribution, and (b) effect of explicit 
spin-dependence on the arrival time distribution of freely falling wavepackets. Section V contains a summary of the 
main results along with certain concluding remarks. 


II. ARRIVAL TIME DISTRIBUTION FOR A TWO-PARTICLE SYSTEM 

Consider a two-body system composed of two non-interacting identical particles in an external field. Identical 
particles are classically distinguishable and obey Maxwell-Boltzmann (MB) statistics, while they are indistinguishable 
in quantum mechanics and obey different statistics. For Fermi-Dirac (Bose-Einstein) statistics the total wavefunction 
must be antisymmetrized (symmetrized) under the exchange of particles in the system. Since particles do not interact, 
solutions of the Schrodinger equation are constructed from two single-particle wavefunctions ipa and pJb as follows @ : 

= 1 pa{zi,t)'lpb{z 2 ,t) , (1) 

'^±{zi,Z 2 ,t) = N±[lpa{zi,t)lljb{z 2 ,t) ±lljb{zi,t)lpa{z 2 ,t)] , (2) 

where upper (lower) sign stands for BE (ED) statistics and the normalization constants are given by N± = [2(1 ± 
Two-body wavefunctions ([2|) are solutions of the two-body Schrodinger equation 

+ + (3) 

where the Hamiltonian is clearly symmetric under the exchange of zi and Z 2 ■ 

The one-body probability density for observing a particle at point z regardless of the position of the other particle, 
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is given by 


pi{z,t) = j dzidz 2 S{z-zi) \'i/{zi,Z 2 ,t)f = j dz 2 \^{z,Z 2 ,t)Y 


= \N±? j 


dZ2 


'lpa{z, t)lpb{z2, t) ± ijjbiz, t)'ll)a{z2, t) 


= \N±\^\\^a{z,t)\'^ ±2'^[{'ll)a{t)\i)bit)) i)l{z,t)i)a[z,t)\ + \'lpb{z,t)\^J , (4) 

where (<)|■!/);,(i)) = J dx ^*(a:, i)^/’{,(a;, t) is the overlap integral and the normalization condition for single-particle 
wavefunctions '0a and tpb has been used. From the unitarity property of the time evolution operator, U{t, 0) = 
one sees that the overlap integral is independent of time, i.e., (0a(t)|0b(t)) = (0a(O)|C/00 0)C/(t, O)|0h(O)) = 
(0a(O)|0b(O)). 

The single-particle continuity equation may be derived [21| in the following way. The partial derivative of pi{z,t) 
with respect to the time is calculated from eq. (|3]) by performing integration by parts and using the fact that the 
wavefunction vanishes at infinity, to obtain 


dpi{z,t) 

dt 


h ^ 
m 


dz 


j dZ2 ^*{z,Z2,t) 


d-^{z,Z 2 ,t) 

dz 


( 5 ) 


Now, using the normalization condition for the single-particle wavefunctions, one obtains the continuity equation 

dpiiz.t) d , , 

T + = “• 

for the one-body probability density pi{z,t) (j4|), with the one-body probability current density given by 

= ^|A^±p3|0*^ -f 06^ ± (0a(t)|0b(t)) 06^ ± (0b(t)|0a(t)) 0*^1 ■ 

For distinguishable particles obeying classical MB statistics, each particle has its corresponding continuity equation, 

( 8 ) 


( 6 ) 


( 7 ) 




with 


and 


j(i)(z,t) = —3<{ 0*(z,t) 
m 


p’-^Hz,t) = |0a(z0)P 

dlpa{z,t) 
dz 


P^‘^\z,t) = |0b(z,t)P , 


j^^\z,t) = —3<{ 0J(z,t) 
m 


dilJb{z,t) 

dz 


So, for a two-particle system, the MB probability density and probability current density are given by 
Pmb(^0) = ^ {p^^\z,t) + p^‘^\z,t)^ , jMB{z,t) = i 


(9) 


( 10 ) 


( 11 ) 


From eqs. m and o, one sees that as long as the single-particle wavefunctions have negligible overlap, i.e., 
(0a(t)|0b(t)) ~ 0, there is no need for symmetrization. As a result, one can ignore indistinguishability of particles 
and thus use MB statistics for which motions of particles are independent. Hence, in the latter case a similar analysis 
to the one-body systems discussed earlier [ij, [T^ may be employed. 

The probability current approach for computation of the arrival time distribution of quantum particles not only 
provides an unambiguous definition of arrival time at the quantum mechanical level |15l- ll8|. b ut also furnishes a way 
of obtaining the classical limit of the mean arrival time for massive quantum particles [22j| . Here we employ the 
probability current approach to study the effect of particle statistics on the arrival time distribution of a two-body 
system. In this approach, the arrival time distribution at a detector location z = Z is given by [H, 


n(z,t) = 




ir dt \ji{Z,t)\ 


( 12 ) 
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As a result, one obtains 


TiZ) = 


dt t n(z, t) 


(13) 


for the mean arrival time at the detector location z = Z. It may be noted that the above expression for the mean 
arrival time follows uniquely in a causal interpretation of quantum mechanics provided by the Bohmian model 
[ 2 ^. Note also that eas. lfT^ and (fT^ represent the single-particle arrival time distribution and the mean arrival time, 
respectively, when the effects of particle statistics of the considered two-body system have been taken into account in 
the computation of the single-particle current density through the process of symmetrization or asymmetrization of 
the wavefunction. In the following analysis we will consider the initial single particle wavefunctions to be Gaussians 
given by 


' 0 i( 2 ,O) = (27rcrgJ i/^exp 


ikiZ — 


(Z - Zci) 


21 




i = a.b 


(14) 


and then compute the arrival time distribution of the cases of first free evolution, and then free fall under gravity. In 
eq-dlll), (Toi, Zci and ki represent respectively the rms width, center and kick momentum of the i-th packet. 


III. EFFECT OF STATISTICS ON FREELY PROPAGATING WAVEPACKETS 


In this section we study the effect of symmetrization and asymmetrization of a two-body freely propagating wave- 
function on the arrival time distribution computed through the probability current approach. The initial one-particle 
wavefunctions are taken to be Gaussian packets given by eq. (1141) . The time evolved wave functions are obtained to be 


'^Az.t) = 


1 


(27rs2Ji/4 


exp • 


2ht 




ht 2 \ ^Oi f 7 I 

^ ^ci ) 1 ^ H” 


2ma. 


Oi 




m 




(15) 


where su = (Toi(l + ih.t/2maQ^) denotes the time-evolved spread of the wavepackets. The overlap integral is given by 


exp 


C'Oa + <^0b 


4(fca - kbYal^alf^ + {Zca - Zcb) + 4f(fca - kb){alf,Zca -I- crl^Zcb) 


and thus, the normalization constants read 


AT ^ _ 1 _ 2 cr 0 aO 0 b 

N± = —=il± - 5 - exp 

V 2 I ^Oa + 


4(CT0a + (^Ib) 


4(fca - hYf^la'^lb + {Zca “ Zcbf 


2 (^§a+^ 0 \) 


- 1/2 


(16) 


(17) 


One sees that the magnitude of the overlap integral decreases exponentially with the separation of the initial wavepack¬ 
ets. As the separation grows, the quantum and classical results approach: 


Pi{z,t) ~ ^{\ilJa{z,t)f + \tj;t{z,t)f) = pMB{z,t) , 


The position expectation value with respect to pi{z,t) reads 
(z) = J zpi{z,t) dz 


(18) 


= i Zc\,a + Zc\,b ± \{lpa{t)\lpb{t))\ 


2 f + hcr^b) - M^caCr^a + ^cbogJX 


I 


m{<yoa + ^Ob) 




(19) 


where Zc/a denotes the classical path of the center of the wavepacket tpa, be., Zc/a = Zca + fikat/m. One sees that the 
expectation value of the center of mass operator is just the position expectation value with respect to pi{z,t), i.e.. 


(^cn.)(t) = (^^)(t) = 


dzidz2 ^^ ^^^ |4'(zi,Z2,t)P 


ijydzizi j dZ 2 \^{zi,Z 2 ,t)\^ + J dz 2 Z 2 J dzi|4'(zi,Z2,t)p| = J dz z pi{z,t) . 


( 20 ) 





















5 


N 



FIG. 1: (Color online) Arrival time distribution at the detector location Z = Q for different initial separation of the wavepackets 
for the free evolution case. Black curves stand for BE statistics, the red ones are for the FD statistics and the green ones show 
MB statistics. Here, ka = kb = —2lao, Zcb = Sctq, cto = 10 ^m, m = nin and Uei = 2mTia%l‘h = 3.165 ms. 


We calculate the arrival time distribution at the detector location Z = 0 by substituting ea. (jl5|) in the expression 
for the probability current given by eq. o, and then substituting the obtained expression in eq. (Ha- The values of 
the parameters used for our numerical computations are given by: doa = cofc = do = 10 /rm, Zca = lOdg, Zch = 8do, 

ka = kb = —2laQ, m = mn = 1.67 x 10“^^ kg, and tref = = 3.165 ms. In figure[T]the arrival time distribution 

has been plotted for different values of initial separation of the wavepackets. It is observed that the inclusion of particle 
statistics through symmetrization and antisymmetrization of a two-body wavefunction has a perceptible impact on 
the arrival time distribution. Such an effect of particle statistics on arrival times of freely propagating particles has 
been studied earlier by using the concept of the crossing state and the second quantization formalism [2^. In the 
present work we further compute the effect of initial separation of the overlapping wavepackets, as is displayed in the 
Figure [TJ We fix the center of and change the center of from one plot to the other. Using eq. (fT6l) and eq. (I13L 
for the parameters of Fig. [T]one obtains 


^caf O'q 

(V'alV'b) 

'^BE j tref 

EpD i^ref 

EMB/tref 

10 

0.6065 

2.371 

2.546 

2.427 

11 

0.3247 

2.518 

2.614 

2.561 

12 

0.1353 

2.682 

2.707 

2.694 

13 

0.04394 

2.826 

2.829 

2.822 


One sees that as the separation of the initial wavepackets grows, the quantum and classical results approach as a 
result of decrement of the overlap integral. 
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IV. VIOLATION OF THE WEAK EQUIVALENCE PRINCIPLE OF QUANTUM MECHANICS 


A. Effect of statistics on freely falling wavepackets 


The arrival time of freely falling wavepackets using the probability current approach has been studied earlier Bin, 
where it was noted that the the arrival time distribution acquires a mass-dependence due to wavepacket spread. In 
the previous section we have seen that the effect of particle statistics impacts the arrival time distribution even for 
the case of free evolution. Now, let us study the effect of symmetrization and asymmetrization of the wavefunction 
of a system of two particles falling freely under gravity. Here again we choose the initial single-particle wavefunctions 
as Gaussians given by eq. da, and set the initial velocity of the particles to be zero, i.e., the particles are dropped 
from rest with ka = kb = 0, and accelerate downwards under gravity with g = 10 m/s^. All numerical calculations 
have been done based on the approximate value 5 = 10 m/s^, since using the exact value of g is not important 
for the arguments and reasoning. The time-evolved single-particle wavefunctions in the uniform gravitational field 
V{z) = mgz are given by, 




(27rs2ji/4 1 2ht 


im 


2 ,2 fit 2 \ f , , Si* 

^ -g* -nr + - [z- zci + — + —r 

12 2,tti(Tq^ J Sti y 2 (jQi 


where su = croi(l -t- ihtl2maQ^). The overlap integral is given by, 


and hence, the normalization constants become 


/ 2crgaCrgb 

{Zca - ZebY ' 

/ 2 , 2 ®xp 

. 4(crg„ -f cr2J_ 


AT ^ 2cr0aCT0b 

N± = — <i 1 ± 9 , .^2 exp 


V 2 I "L 2 (u 2 „+Uo\)J 


{Zca Zcb) 


- 1/2 


, ( 21 ) 


( 22 ) 


(23) 


One sees that both overlap integral and normalization constants are independent of g. Here again the magnitude of 
the overlap integral decreases exponentially with the separation of the initial wavepackets, and as in the case of freely 
evolving wavepackets, the quantum results of the probability density and current approach the classical expressions 
as the separation grows. The position expectation value is given by 


(2)± = J zpi{z,t) dz 

= ^ ( 34 ) 

L V ^ ^Oa + '^Ob / J 

where Zci,a shows the classical path of the center of the wavepacket ipa, be, n,a = Zca — gt^ The position expectation 
value with respect to pi{z, t) is the expectation value of the center of mass operator. The quantity {z)± is independent 
of mass. 

For the free fall case and for croa = crob = one obtains 



(25) 


for the width of the single-particle distribution. One clearly sees this quantity is independent of the gravity strength 
g and its growth with time is slow for large masses. 

We compute the arrival time distribution at the detector location Z = 0 by substituting the expression for the time- 
evolved wavefunction ea. (| 2 T|) in the expression for the probability current given by eq.([7]), and then using ea. (|T^ . 
Using the three quantities fi, g and uo, one can construct a quantity mg = fi/^/g^ with the dimension of mass, 
which has the value toq = 1.055 x 10“^^kg k, m„, the neutron mass. We use to make the mass m dimensionless. 
The values of the parameters used for our numerical calculations are given by: coa = <^ob = erg = 10 /rm, Zca = lOcrg, 
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FIG. 2: (Color online) Arrival time distribution at the detector location Z = Q for different initial separation of the wavepackets 
for the free fall case. Black curves show BE statistics, the red ones are for the FD statistics and the green ones show MB 
statistics. Here, z^h = Sao, ao = 10 /rm, g — 10 m/s^, m = m„ and trei ~ 2m^rj%l‘h — 3.165 ms. 


Zch = StToj ^ = 1-67 X 10“^^ kg, and fret = = 3.165ms. In figure [5] the arrival time distribution has 

been plotted for different values of initial separation of the wavepackets. In this case too the particle statistics has an 
impact on the arrival time distribution. Using eq. (12211 and eq. da, for the parameters of Fig. [2] one obtains 


^ca! 

(V'alV'fe) 

TBE/iref 

Tpo/tref 

TMB/tref 

10 

0.6065 

1.339 

1.341 

1.340 

11 

0.3247 

1.374 

1.375 

1.374 

12 

0.1353 

1.407 

1.407 

1.407 

13 

0.04394 

1.439 

1.439 

1.439 


Such an effect of symmetrization and asymmetrization of a two-body wavefunction on arrival times of freely falling 
wavepackets may be regarded as nonlocal (in the sense that the single-particle arrival time distribution depends on 
the spatially separated second particle, as well), and thus contrary to the tenet of the local weak equivalence principle 
of classical gravity, which forms the inspiration of the statement of WEQ. 

Now, in order to show the explicit mass dependence of arrival time for particles in free fall we plot the mean arrival 
time versus mass in the figure [3l Here again, the mean arrival time is computed at the detector location Z = 0 using 
eq. (113 and eg. (1131) after substituting the expression for the time-evolved wavefunction eg. (1211) in the expression for 
the probability current given by eq. (CD- The values of the parameters used are as before. Curves with label ilja {'4’b) 
show mean arrival time for wavepackets ipa {'4’b) respectively. The MB mean arrival time is given by ttvib = . 

One sees that for all types of statistics the mean arrival time decreases with mass at first and then becomes constant 
for large mass. The effect of symmetrization and asymmetrization of the wavefunction has a perceptible effect on 
the magnitude of the violation of WEQ. However, the violation of WEQ observed explicitly for low masses tends to 
disappear smoothly in the limit of large mass. Hence, even with the inclusion of particle statistics, the classical limit 
emerges smoothly through this approach. 

One sees each of the three statistics gives a different mean arrival time in the large mass limit. The difference 
between the mean arrival times of bosons and fermions is related to the difference between the two quantities (i) 
position expectation value, with respect to the one-body probability density pi{z, t) and (ii) spread in position. 




















FIG. 3: (Color online) Mean arrival time at the detector location Z = Q versus mass for the free fall case. Here, ka = kh = 0, 
Zca = 10(70, Zcb = 8(70, (70 = 10 fJ,m and tref = 2mn(7o/?i = 3.165 ms. In b) we have included mean arrival time of single-particle 
systems which are described by wavepackets and t/jb, i.e., Ta and Tb- 


Az(t) = •\/(i2y3'(^y2^ At time t = t±, 


t± 


^ca H” ^cb 

2 =F exp 

{Zca-ZcbY 

4(Tg 

9 

2 ± exp 

(^ca ■2cb)^ 

4a^ 


(26) 


the center of the one-body probability density pi(z,0) arrives at the ground, the detector location. We have plotted 
(^)±(0) and Az±{t±) as a function of mass in Fig. H) As this figure shows Az+{t+) and Az-{t-) corresponding to 
bosonic and fermionic statistics, respectively, do not approach each other even in the large mass limit. 


B. Spin-dependent arrival time distribution 

In this subsection we discuss the effect of explicit inclusion of spin in the probability current density on the arrival 
time distribution of particles in free fall. It has been noted earlier that in nonrelativistic quantum mechanics the 
form of the probability current density is not unique. It was shown [l^ that the Dirac equation implies a unique 
expression for the probability current for spin-1/2 particles in the nonrelativistic limit. Uniqueness of the probability 
current is a generic consequence of any relativistic quantum dynamics. Effects of spin on the arrival time distribution 
for freely evolving wavepackets have been computed for fermions [l^ as well as for bosons [l^. Here we obtain the 
spin-dependent arrival time distribution for a symmetric two-dimensional Gaussian wavepacket in order to compare the 
effect of explicit presence of spin in the probability current density with the effect of particle statistics for wavepackets 
falling freely under gravity. Taking into account the spin contribution to the probability current density for particles 
with spin-1/2, one obtains (l^ 


j(x,t;s) 


jsch(x,t) -k — X s 

m 

t)V'0(x, <)]-I-t)V'0(x, t)] X §} , 
m 


(27) 

(28) 
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FIG. 4: (Color online) (2)±(0) and Az±{t±) for the free fall case. Black curve shows BE statistics and the red one is for the 
FD statistics. Here, z^a = lOao z^b = Scto, fro = 10 /rm and g = 10 m/s^. 


for a spin eigenstate x in the absence of a magnetic field. Here, jsch(x, t) is the usual Schrodinger current, whereas 
the second term denotes the spin-dependent contribution with s = For ^ two-dimensional system with a 

wavefunction which is in the factorized form ^(x, t) = ' 0 a;(a:, t)' 02 (z, t), the probability current density takes the form 


m 

m 


(29) 

(30) 


for s = (0,1,0). In the uniform gravitational field V{z) = mgz, by taking the initial wavefunction as a two-dimensional 
Gaussian wavepacket, 

(z - z,)2 


’fl^oix.z) = 




4cr2 


■ ikoz 


the modulus of the Schrodinger current takes the form 

+ f3{t)m^ + f 2 {z, t)m‘^ + fi{z, t)m + /o(x, 2 , t) 


|jsch(^5 -^5 I - 


IGnm'^aQaf 


X exp < — 


+ [z - Zcl(t)]^ 

2af 


while the modulus of the total spin-dependent probability current density is given by 

X^ + [z- Zcl(t)]^ 


, .x, _ + hi{x, t)m + ho{x, z, t) 

SJI o_™_ ..^3 OXP) 20-2 


Snmaoaf 


where. 


foix, z, i) = + 4[a:^ + {z - z^)^] -h ^gt^{-z + z^)} , 

= -lQn^koaQt{gt^-2z + 2zc) , 

/ 2 (z,t) = lQh^(jQ[Akla^+gt^{gt^-2z + 2zc)\ , 
fsit) = - 128 ghkoCT^t , 
fi{t) = Mg'^alt^ , 


(31) 


(32) 


(33) 


and 


ho{x,z,t) = h^llOklaQ + g'^t'^ - lOkoa^x + A[x'^ + {z - Zc)^]+Agt'^{-z + Zc)} , 
hi{x,t) = l(Sg%aQt{—2kQcf^ + x) , 
h2{t) = le^^CTpt^ . 
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FIG. 5: (Color online) a) Spin-dependent (black points) and Spin-independent (red points) mean arrival time at the detector 
plane Z = 0 and b) their difference, versus mass for the free fall of the 2D symmetric Gaussian packet (I31II . The parameters of 
the falling wavepacket have been chosen as Zc = Sag and ctq = 10 ^tm, rUn = 1.67 x 10~^^ kg and tret = 2mnO'o/^ = 3.161 ms. 


In eqs. (1321) and (1551) . at = goy^l + and Zc\{t) = —^gt^ + Uzt + Zc with Uz = hko/m. 

In order to obtain the arrival time distribution at the plane z = Z we integrate over the space coordinate x. Thus, 
the normalized arrival time distribution is given by 


nsch(^, t) 


|jsch(a;,^,0l 

\isch{x, Z,t)\ 


in the absence of spin contribution, while it reads 


U{Z,t;s) 


ir dtJZodx \i{x,Z,t;s)\ ’ 


(34) 


(35) 


considering the spin effect. We have integrated over x from — oo to oo to compare the results with the corresponding 
one-dimensional system. The corresponding mean arrival times are respectively given by 


TSchC^) 

t(Z;s) 


poo 

/ dttILsch{Z,t), 

Jo 

pOO 

/ dt t n(Z, t; s) . 
Jo 


(36) 

(37) 


In the Figure [SI the mean arrival times rsch and r(s) at detector location Z = 0 have been plotted versus mass for 
free fall from rest ( fcp = 0). The parameters of the initial wavepacket have been chosen as Zc = SgQ and gQ = 10 /rm, 

TOn = 1.67 X 10“^^ kg and tref = = 3.161 ms. One sees that the mean arrival time is mass-dependent, thus 

signifying a violation of WEQ. Inclusion of spin causes the mean arrival time to decrease with a small magnitude 
compared to the case of no spin. Hence, in comparison with the results for the case considering the effect of particle 
statistics without spin-dependence displayed in the figure |3l we find that the contribution of including spin explicitly 
at the level of the probability current leads to a much smaller modification of the mean arrival time values. However, 
here (Figure jS]) again we find that the mass dependence of the mean arrival time vanishes for large mass, signifying 
emergence of WEQ in the classical limit. 
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V. CONCLUSIONS 

In this work we have investigated the effect of particle statistics on the arrival time distribution of Gaussian wave 
packets. We have employed the probability current approach [Tsi l for computing the arrival time distribution for a 
system of two particles. The probability current approach provides an unambiguous definition [Tsl of the arrival 
time of quantum particles, and also leads to a proper classical limit of quantum dynamics. The single-particle 
arrival time distributions using the probability current approach have been studied earlier in the literature for the 
cases of free evolution, leading to the interesting effect of the spin-dependence [3, [iM of arrival time. It was also 
observed earlier that the mean arrival time for particles freely falling under gravity acquires mass-dependence, thus 
signifying violation of the weak equivalence principle of quantum mechanics |12l . [l3| . Such violation stems essentially 
from the spread of propagating wavepackets that is mass-dependent. 

In the present analysis we have shown for the first time that symmetrization or asymmetrization of a two body 
wavefunction leads to modification of the arrival time distribution even for freely propagating wavepackets. The impact 
of particle statistics on the arrival time distribution is observed for initially overlapping Gaussian wavepackets. For 
the case of particles falling freely under gravity, the mean arrival time acquires mass dependence signifying violation 
of WEQ. The magnitude of violation is seen to be different for MB, BE and ED statistics as a consequence of the 
difference in the mass dependence of the mean arrival time for the case of different statistics. We have also compared 
such a violation under the influence of particle statistics with the violation of WEQ obtained through the explicit 
introduction of particle spin in the probability current. The modification of the mean arrival time in the former 
approach is shown to be larger than in the latter for the case of Gaussian wavepackets considered here. However, 
in both the approaches, WEQ emerges smoothly in the limit of large mass. It may be noted that each of the three 
statistics gives a different mean arrival time in the large mass limit. Such difference stems essentially from the 
difference in the position expectation value and the spread of the wave packets for the respective statistics, and may 
be viewed as a remnant effect of particle statistics in the classical limit. 

We conclude by observing that the effect of particle statistics in different quantum phenomena is regarded to be of 
much importance motivating newer tests beyond traditional arenas [^ . Our approach of studying the influence of 
particle statistics on the arrival time distribution provides an independent manifestation of the fundamental quantum 
property of indistinguishability. Moreover, the analysis of transit and flight times iEl is a key ingredient in 
experiments involving Bose condensates which form a major avenue of revealing the effect of quantum statistics. 
The results reported in the present paper using symmetric Gaussian wavepackets furnish an example of the effect of 
particle statistics on the arrival time distribution. Though the magnitude of the difference in observable quantities 
resulting from different statistics, as computed in the present work may not be high enough for direct experimental 
detection, they provide an in principle signature of differential violation of WEQ for different statistics. More elaborate 
calculations are needed using different sets of parameters, and various forms of wavepackets such as asymmetric fisj ] 
or non-Gaussian [T^ ones in order to raise the possibility of observational verification of such effects. 
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